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Abstract 


Wigley  hulls  have  been  widely  used  as  test  cases  for  hydrodynamic  computer  programs  because 
they  can  be  modelled  exactly  in  the  computer,  and  because  there  are  extensive  experimental 
data  describing  the  flow  around  them.  In  this  communication,  it  is  shown  how  Wigley  hulls  can 
be  described  using  the  HLLSRF  hull  representation  system.  A  user’s  guide  for  a  Fortran  77 
program  which  generates  a  HLLSRF  data  file  for  a  Wigley  hull  is  included. 

L A,  > -x )  <jjyLr^ 


/r. 


i- 


Resumd 


(  K  P- 


Les  coques  Wigley  ont  6t6  largement  utilises  pour  l'essai  des  programmes 
hydrodynamiques  informatis£s,  car  elles  peuvent  etre  modeiis£es  avec  exactitude  par 
l'ordinateur  et  parce  qu'on  dispose  d’une  grande  masse  de  donndes  exp^rimentales  relatives 
ci  l'dcoulement  d'eau  autour  d'elles.  Dans  la  pr£sente  communication,  on  montre  comment 
les  coques  Wigley  peuvent  etre  decrites  au  moyen  du  systfcme  HLLSRF  de  representation 
des  coques.  Est  inclus  un  guide  d'utilisation  du  programme  Fortran  77  qui  genfcre  les 
donndes  d'une  coque  Wigley  dans  le  systfcme  HLLSRF. 
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Notation 


a  hull  form  parameter 
B  the  hull  half-breadth 

B^\  the  nth  B-spline  of  order  kx  with  respect  to  the  knot  sequence 

B(n%  the  nth  B-spline  of  order  ky  with  respect  to  the  knot  sequence  t ^ 
D  the  hull  depth 
Cb  block  coefficient 
kx  the  order  of  the  B-splines  B^ki 

ky  the  order  of  the  B-splines 
L  length  between  perpendiculars 
Nx  the  number  of  B-splines  B^l 
Ny  the  number  of  B-splines  B^ 
p  —TjD 

t the  knot  sequence  used  to  define  B^kx 

t[y^  the  knot  sequence  used  to  define  B^ky 
T  the  hull  draft 

x  HLLSRF  non-cartesian  coordinate 
xap  the  station  number  of  the  aft  perpendicular 
xpp  the  station  number  of  the  forward  perpendicular 
X  HLLSRF  cartesian  coordinate 
X'  cartesian  coordinate  used  to  define  the  Wigley  hull 
y  HLLSRF  non-cartesian  coordinate 
Y  HLLSRF  cartesian  coordinate 
Y'  cartesian  coordinate  used  to  define  the  Wigley  hull 
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Z  HLLSRF  cartesian  coordinate 
Z'  cartesian  coordinate  used  to  define  the  Wigley  hull 
Qjn  spline  coefficients  used  to  define  Y 
{3jn  spline  coefficients  used  to  define  Z 
7„  coefficients  defined  in  equation  (2.4) 

A„  coefficients  defined  in  equation  (2.4) 

Hn  coefficients  defined  in  equation  (2.4) 
rn  coefficients  defined  in  equation  (2.4) 
fy(x)  function  defined  in  equation  (3.9) 

<7y(x)  function  defined  in  equation  (3.10) 
fz{x)  function  defined  in  equation  (3.11) 
gz(x)  function  defined  in  equation  (3.12) 


vi 


1  Introduction 


Because  Wigley  hulls  can  be  represented  exactly  by  computers,  they  have  been  widely 
used  as  test  cases  of  hydrodynamic  programs,  particularly  for  calculations  of  wave  resistance 
and  boundary  layer  flow.  Experimental  data  on  the  flow  about  Wigley  hulls  are  extensive. 

The  hull  representation  system  HLLSRF1 ,  developed  at  DREA  for  use  in  hydrodynamic 
computations,  is  capable  of  representing  Wigley  hulls  exactly;  however,  the  correct  form  for 
that  representation  is  not  readily  apparent.  In  this  communication  the  HLLSRF  representation 
of  the  Wigley  hull  is  derived. 

Appendix  A  is  a  user’s  guide  for  a  Fortran  77  program,  WIGLEY,  which  implements  the 
calculations  described  here. 


2  The  HLLSRF  Hull  Representation 


HLLSRF  uses  B-splines  of  arbitrary  order  to  represent  hulls  (see  de  Boor2  for  further 
information  on  B-splines).  The  hull  representation  defines  a  functional  form  for  the  cartesian 
coordinates  on  the  hull  surface  and  two  non-cartesian  hull  coordinates,  x  and  y.  The  cartesian 
coordinates  are  denoted  X,  Y,  and  Z;  X  varies  between  0  and  L  as  one  goes  from  bow  to  stem; 
Y  is  the  distance  from  the  centreplane;  and  Z  is  the  distance  above  the  baseline.  The  HLLSRF 
representation  relates  (A',Y,  Z)  to  (x,y)  as  follows. 


X(x,y) 

Y{*,y) 


Z(x,y) 


( x  -  xp-p)L 
zap  -  xfp 


j= 1 n=l 


j=l n=l 


(2.1) 

(2.2) 

(2.3) 


where 

L  is  the  length  between  perpendiculars, 

B  is  the  half-breadth, 

D  is  the  hull  depth, 

xpp  is  the  station  number  at  the  forward  perpendicular, 
xap  is  the  station  number  at  the  aft  perpendicular, 
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Bn,lx(x)  is  a  B-spline  of  order  kx  corresponding  to  a  knot  sequence  t\x\i  =  1, . , 
each  n  —  1, . . . , Nx, 

. . ,  Nx  +  kx,  for 

Bi%(y)  is  a  B-spline  of  order  ky  corresponding  to  a  knot  sequence  t\v\i  =  1, . , 
each  n  =  1,.. . , Ny,  and 

..,Ny  +  ky,  for 

ajn,0jn  are  coefficients. 

If  the  spline  coefficients  ojn  and  f3jn  have  the  form 

ajn  =  7 j^n  i  &jn  —  Tn 

(2.4) 

then  the  HLLSRF  representations  of  the  Y  and  Z  coordinates  have  the  form 

Y  =  BfY(x)gY(y ) 

Z  =  Dfz{x)gz{y) 

(2.5) 

(2-6; 

where 

fY(X)  =  y:kb[xI(x) 

(2.7) 

Nv 

fjy(y)  =  Yl'riBjy,)kv(y ) 

(2.8) 

H 

SI- 

03 

.5 

I 

II 

(2.9) 

Ny 

9z{y)  = 

j=i 

(2.10) 

HLLSRF  can  represent  exactly  any  hull  which  can  be  written  in  the  form  of  equations  (2.1), 
(2.5),  and  (2.6)  where  fY,  fz,  9Y »  and  gz  are  piecewise  polynomials. 


3  Representation  of  Standard  Wigley  Hulls 


The  class  of  mathematical  hull  forms  known  as  Wigley  hulls  is  usually  defined  as  follows3. 

v  (,  z'2\  ( ,  4*,2\  A  4x,2\  , 

B  ~  \l  T2 )  L2  )  (1  +  a  L2  )  (3-1) 

Here  T  is  the  draft  and  X ',  Y'  and  Z'  are  cartesian  coordinates  such  that  X'  varies  between 
—  Lf 2  and  L/2  as  one  goes  from  bow  to  stem;  Y'  is  the  distance  from  the  centreplane;  and  Z' 
is  the  distance  below  the  waterline.  In  terms  of  the  HLLSRF  cartesian  coordinates 


X' 

Y' 

Z' 


X-L/2 

(3-2) 

Y 

(3.3) 

D-Z 

(3-4) 

2 


Setting  D  =  T,  equation  (3.1)  can  be  rewritten 


Y 

B 


z  ( 

1 4X  ! 

,  x\ 

1  ,  ( 

1  2*V1 

T  V 

0 

t) 

‘-t) 

1  +  0,  ( 

l-T)\ 

(3.5) 


The  hull  form  parameter  a  may  be  used  to  vary  the  block  coefficient,  Cb,  which  is  related 
to  a  by 

c»  -  5  0 + !)  (3-6) 

The  value  of  a  must  be  in  the  range  (  —  1,1)  to  produce  realistic  hull  forms;  for  a  <  —1,  the 
value  of  Y  will  be  negative  at  some  stations,  while  if  a  >  1,  the  value  of  Y  will  exceed  B  at  some 
stations.  Hence,  the  block  coefficient  must  be  in  the  range  (16/45,24/45)  =  (0.35555,0.53333). 


Let  xpp  =  — 1,  xap  =  1,  and  y  =  Zj D.  Then  from  equation  (2.1)  one  has  x  —  2X/ L  —  1  = 
X' / L  and  equation  (3.5)  becomes 

Y_ 

B 

£ 

D 


=  y(2  -  J/)(l  -  2:2)(1  +  ax2) 

(3.7) 

=  y 

(3.8) 

ns  (2.5)  and  (2.6)  with 

/y(x)  =  (1  -  x2)(l  +  ax2) 

(3.9) 

9y(v)  =  y(2-y) 

(3.10) 

fz{x)  =  1 

(3.11) 

gz(y)  =  y 

(3.12) 

Both  fy  and  fz  are  polynomials  of  fifth  order  or  lower  over  the  region  [—1,1];  therefore, 
they  can  be  represented  by  B-splines  generated  by  setting  Nx  =  5,  kx  =  5,  and  the  x-knot 


sequence,  t\x\  to  be  {  — 1,  — 1,  —  1,  —  1,  —  1, 1, 1 

BiM  = 
£2,5(1)  = 

£3,5(2;)  = 
£4,5(1)  = 
£5,5(2;)  = 


1,1,1}.  The  B-splines  in  x  are 

(1  -  x)4 

16 

(3.13) 

(1  +  i)(1-i)3 

4 

(3.14) 

3(1  -x2)2 

8 

(3.15) 

(1-xXl+i)3 

4 

(3.16) 

(1  +  x)4 

16 

(3.17) 

(See  de  Boor2,  chapter  IX,  for  methods  of  generating  B-splines  from  their  associated  knot 
sequences.)  Equations  (2.7),  (2.9),  (3.9)  and  (3.11)  are  satisfied  by  setting 


Aj  =  A5  =  0  ;  Aj  =  A4  =  1  +  a 


A3  = 


4(1  ~a) 

3 


(3.18) 


3 


t\  =  t2  =  r3  =  r4  =  r5  =  1 


(3.19) 


Both  gy  and  gz  are  polynomials  of  'bird  order  or  lower  over  the  region  [0,1].  Moreover, 
each  vanishes  at  y  =  0,  so  that  they  can  each  be  represented  by  setting  Ny  =  2,  kx  =  3,  and 
the  y-knot  sequence,  t\y\  to  be  {0,0,1, 1,1}.  The  B-splines  in  y  are 

£1,3(1/)  =  2y(l  -  y)  (3.20) 

B2,3(y)  =  y2  (3.21) 

Equations  (2.8),  (2.10),  (3.10)  and  (3.12)  are  satisfied  by  setting 

7i  =  72  =  1  (3.22) 

^  =  1/2  ;  n2  =  1  (3.23) 

The  corresponding  values  of  atjn  and  /?jn  are  given  in  Table  1. 


®jn 

n  =  1 

n  =  2 

n  =  3 

n  =  4 

n  =  5 

j  =  1 

0 

1  "f  fl 

4(1  -  a)/3 

1  +  a 

0 

3=  2 

0 

1  -J-  CL 

4(1  -  a)/3 

1  *4"  ci 

0 

0jn 

n  =  1 

n  =  2 

n  =  3 

n  —  4 

n  =  5 

j  =  1 

1/2 

1/2 

1/2 

1/2 

1/2 

3=  2 

1 

1 

1 

1 

1 

Table  1:  ajn  and  0jn  for  Wigley  hulls  below  the  waterline  only. 


4  Wigley  Hulls  with  Freeboard 


For  hydrodynamic  calculations  where  the  water  surface  is  free,  it  is  often  useful  to  extend 
the  standard  Wigley  hull  above  the  waterline.  This  may  be  done  by  extending  each  station 
vertically.  Figure  1  shows  an  offset  diagram  of  a  Wigley  hull  extended  in  this  fashion;  the 
additional  horizontal  line  is  the  waterline. 


The  formula  for  the  extended  hull  corresponding  to  equation  (3.5)  is 


y 

B 


?K)h(-¥)1 


if  Z  <  T 
if  Z  >  T 


(4.1) 


4 


0.13  L 


0.15  L 


Figure  1:  Offset  diagram  of  a  Wigley  hull  with  freeboard 


Retaining  the  identification  y  =  Z/D  so  that  y  varies  between  0  and  1,  the  counterpart  of 
equation  (3.7)  is 

V  f  -  (2-  -)  (1  -  z2)(l  +  ax2)  if  g  <  p 

g  =  \  V'  (4.2) 

(  (1  -  X2)(l  +  ax2)  if  y>p 

where  p  =  T/D.  This  equation  is  again  of  the  form  of  equations  (2.5)  and  (2.6)  with  fy  and 
fx  unchanged  from  the  case  with  no  freeboard.  The  functions  gy  and  gz  are  given  by 


9r(y) 

9z{y) 


(  (1  -  i2)(l  +  ax2)  if  y  >  p 

y 


(4.3) 

(4.4) 


Both  gy  and  gz  are  piecewise  polynomials  of  order  three  or  less  over  the  region  [0, 1]  with  a 
single  breakpoint  at  y  =  p.  Since  each  vanishes  at  y  =  0,  they  can  be  represented  by  setting 
Ny  =  3,  kx  =  3,  and  the  y-knot  sequence,  t\v\  to  be  {0,0, p,  1, 1, 1}.  The  B-splines  in  y  become 


#1,3  (y) 


f  ^-(l  +  rtv)  fojy<p 

I  p* 

(±^aZ  {o,s>p 

\  1  -p 


(4.5) 
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BiM  = 


B3M  =  { 


r 

p 


for  y  <  p 


(l-y)((2  -p)y-p)  , 

— (wji —  «»t>p 
0  for  y  <  p 

for  y  >  p 


l  (1-P)2 

Equations  (2.8),  (2.10),  (4.3)  and  (4.4)  are  satisfied  by  setting 


(4.6) 


(4.7) 


7i  =  72  =  73  =  1  (4.8) 

Pi=p/ 2  ;  H2  =  (l  +  p)/2  ;  M3  =  l  (4.9) 

The  corresponding  values  of  ajn  and  (3jn  for  the  Wigley  hull  with  freeboard  are  given  in  Table  2. 


n  =  1 

n  =  2 

n  =  3 

n  =  4 

n  =  5 

3  =  1 

0 

1  -f*  d 

4(1  —  a)/3 

1  +  a 

0 

IB1 

0 

1  +  a 

4(1  —  a)/3 

1  +  CL 

0 

IEBI 

0 

1  +  a 

4(1  —  a)/3 

1  +  0 

0 

0jn 

n  =  1 

n  =  2 

n  =  3 

n  =  4 

n  =  5 

3  =  1 

p/2 

p/2 

p/2 

p/2 

p/2 

033 

(l  +  P)/2 

(l  +  P)/2 

(l  +  P)/2 

(1  +  P)/2 

(1  +  P)/2 

IEBI 

1 

1 

1 

1 

1 

Table  2:  ajn  and  (3jn  for  Wigley  hulls  with  freeboard. 


5  The  Quadratic  Wigley  Hull 

The  Wigley  hull  for  which  a  =  0  is  especially  simple  because  Y  is  quadratic,  not  quartic, 
in  X'/L.  It  is  sometimes  called  the  quadratic  Wigley  hull.  Because  Y  is  quadratic  in  X' / L 
and  hence  in  x,  the  i-B-splines  need  only  be  of  third  order  rather  than  fifth  order.  Because 
B-splines  become  more  costly  to  evaluate  as  the  order  increases,  and  because  this  particular 
form  of  the  Wigley  hull  is  widely  used,  it  is  worth  generating  a  third  order  representation  of  it. 
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Since  the  dependence  of  Y  upon  y  is  unchanged,  it  is  only  necessary  to  consider  a  modi¬ 
fication  in  the  B-splines  which  depend  on  x.  From  equations  (3.9)  and  (3.11)  with  a  =  0,  one 
has 


fY(x)  =  (1  -  X 2) 

(5.1) 

fz{x)  =  1 

(5.2) 

(5.3) 

These  functions  must  be  represented  in  terms  of  B-splines  of  order  three  in  x. 
kx  =  3,  and  the  knot  sequence  t\x^  to  be  {  — 1, -1,  —  1, 1, 1, 1),  one  has 

Setting  Nx  =  3, 

- 

BM  =  ^ 

(5.4) 

«.)  =  (1v2) 

(5.5) 

=  (ir)2 

(5.6) 

Equations  2.7,  (2.9),  (5.1)  and  (5.2)  are  satisfied  by  setting 

Ax  =  A3  =  0  ;  A2  =  2  (5.7) 

n  =  r2  =  r3  =  1  (5.8) 

The  corresponding  values  of  ajn  and  (3jn  for  the  Wigley  hull  without  freeboard  and  with 
freeboard  are  given  in  Tables  3  and  4. 


6  Concluding  Remarks 

The  spline  coefficients  necessary  to  represent  Wigley  hulls  using  HLLSRF  have  been  de¬ 
rived.  The  Wigley  hulls  can  be  represented  with  or  without  freeboard.  Once  Wigley  hulls  have 
been  represented  in  this  form,  the  HLLFLO4  suite  of  hydrodynamic  programs  can  be  used  to 
predict  the  flow  around  them. 


4 
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Table  3:  ajn  and  for  quadratic  Wigley  hulls  below  the  waterline  only. 


Table  4:  ajn  and  for  quadratic  Wigley  hulls  with  freeboard. 
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Appendix  A  User’s  Guide  for  Program  WIGLEY 


Purpose:  The  program  WIGLEY  creates  a  HLLSRF  data  file  for  a  Wigley  hull  with  or 
without  freeboard. 


Input:  All  input  to  the  program  is  entered  interactively  by  the  user.  When  the  program 
starts,  the  user  is  prompted  for  the  hull  form  parameter,  a. 

Enter  the  hull  form  parameter  => 

The  hull  form  parameter  must  be  in  the  range  (-1.0, 1.0).  If  it  is  not,  an  error  message  will 
be  written  and  the  user  reprompted.  When  a  valid  hull  form  parameter  has  been  entered,  the 
corresponding  block  coefficient  (see  equation  (3.6))  is  reported  to  the  user. 

Next  the  user  is  prompted  for  the  units  in  which  the  hull  dimensions  are  measured. 
Enter  the  units  for  the  hull  dimensions  ■> 

Any  string  of  ten  or  fewer  characters  will  do:  e.g.  ‘ft’,  ‘metres’,  or  ‘  ’. 

Next  the  hull  dimensions  themselves  are  required. 

Enter  the  length  between  perpendiculars,  the  half -breadth,  the  depth, 
and  the  draft  => 

Each  of  the  hull  dimensions  must  be  positive  and  the  draft  must  not  exceed  the  depth.  If  there 
is  an  error  in  input,  an  error  message  will  be  sent  and  the  user  reprompted. 

These  are  all  the  data  necessary  to  define  the  Wigley  hull. 

Output:  WIGLEY  creates  an  ASCII  file  WIGLEY.HLL  which  contains  the  HLLSRF  repre¬ 
sentation  of  the  hull. 

Usage:  On  the  DREA  DEC-20/60,  WIGLEY  can  be  run  by  typing  the  following  commands. 

DECLARE  PCL  PS : <HALLY>HLLFL0 , PS : <HALLY>WIGLEY 
WIGLEY 

The  first  command  causes  the  command  WIGLEY  to  be  defined  in  the  EXEC.  Typing  WIGLEY 
then  causes  the  program  to  be  run.  The  program  can  be  run  again  simply  by  repeating  the 
command  WIGLEY. 
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